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Abstract 

The non-exponential relaxation is shown to result from subordination by inverse 
tempered a-stable processes. The main feature of tempered a-stable processes is 
a finiteness of their moments, and the class of random processes includes ordinary 
a-stable processes as a particular case. Starting with the Cole-Cole response this 
subordination approach establishes its direct link with the Cole-Davidson law. We 
derive the relaxation function describing the tempered relaxation. The meaning of 
the empirical response function is clarified. 
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1 Introduction 



The major feature of dynamical processes in many complex relaxing systems 
is their stochastic background [T|2] . Particularly, in any dielectric (complex) 
system under an week external electric field (external action) only a part (ac- 
tive dipoles or objects) of the total number of dipoles is directly governed by 
changes of the field. But even those dipoles, not contributing to the relax- 
ation dynamics, can have an effect on the behavior of active dipoles [3]. If the 
dipoles interact with each other, then their evolution has a random character. 
Consequently, the behavior of such a relaxing system as a whole will not be 
exponential in nature. In this case the macroscopic behavior of the complex 
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systems is governed by "averaging principles" like the law of large numbers 
following from the theory of probability [4]. The macroscopic dynamics of 
complex systems is not attributed to any particular object taken from those 
forming the systems. The finding out an "averaged" object representity for the 
entire relaxing system is not simple. The relation between the local random 
characteristics of complex systems and the universal deterministic empirical 
laws requires a probabilistic justification. There are some points of view on 
this problem. One of well-known them is based on randomizing the parame- 
ters of distributions that describes the relaxation rates in disordered systems. 
With regard to the dielectric relaxation, each individual dipole in a complex 
system relaxes exponentially, but their relaxation rates are different and obey 
a probability distribution (continuous function) [3f5] . This approach is succes- 
sive for getting many empirical response laws and their classification, but it 
sometime becomes enough complicated to interpret their interrelations and to 
derive macroscopic response equations.. 

In this paper we suggest an alternative approach to the analysis of non- 
exponential relaxation. It is based on subordination of random processes. Re- 
call that in the theory of anomalous diffusion the notation of subordination 
occupies one of the most important places (see, for example, [6] and references 
therein). So, a subordinated process Y(U(t)) is obtained by randomizing the 
time clock of a random (parent) process Y(t) by means of a random process 
U (t) called the directing process. The latter process is also often referred to as 
the randomized time or operational time [7|. In the common case the process 
Y may be both random and deterministic in nature. The subordination of 
random processes is a starting point for the anomalous diffusion theory. 

We develop this approach to relaxation processes. It gives an efficient method 
for calculating the dynamical evaluating averages of the relaxation processes. 
In this connection Section [2] is devoted a presentation of recent achievements 
of this method. Starting with the description of the two-state system evolution 
as a Markovian process, we develop the analysis on subordinated random pro- 
cesses. The processes differ from the Markovian ones by the temporal variable 
becoming random. In this context the Cole- Cole relaxation is an evident exam- 
ple. In Section [3] we consider the tempered a-stable processes. They overcome 
the infinite-moment difficulty of the usual (not tempered) a-stable processes. 
As applied to the anomalous diffusion, the tempering gives a preserving the 
subdiffusive behavior for short times whereas for long times the diffusion is 
something like normal. Using the processes in Section m we develop a sub- 
ordination scheme for the description of the tempered relaxation. Section [5] 
formulates major properties of such relaxation. We show that it has a direct 
relation to the well-known experimental laws of relaxation, in particular, to 
the Cole-Davidson law. Finally, the conclusions are drawn in Section 
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2 Relaxation in Two-State Systems 

The simplest ordinary interpretation of relaxation phenomena is based on the 
concept of a system of independent exponentially relaxing objects (for exam- 
ple, dipoles) with different (independent) relaxation rates [8]. The relaxation 
process, following this law (called Debye's), may be represented by behavior 
of a two-state system. Let N be the common number of dipoles in a dielectric 
system. If Nf is the number of dipoles in the state f, N± is the number of 
dipoles in the state 4- so that N = + N±. Assume that for t = the system 
is stated in order so that the states t dominate, namely 

NM = 0) , s Ni(t = 0) 

n N , =n r {0) = l, l \ ; =n ; (0) = 0, 

where is the part of dipoles in the state T> n± the part in the state I. Denote 
the transition rate by w defined from microscopic properties of the system (for 
instance, according to the given Hamiltonian of interaction and the Fermi's 
golden rule). In the simplest case (D relaxation) the kinetic equation takes the 
form 

h t (t) -w{n±{t) -n t (t)} = 0, 

h l (t)-w{n t (t)-n i (t)} = 0, 1 1 

where, as usual, the dotted symbol means the first-order derivative. The re- 
laxation function for the two-state system is 

D (t) = 1 - 2n±{t) = 2n t (t) - 1 = exp(-2wt). 

It is easy see that the steady state of the system corresponds to equilibrium 
with n^(oo) = n^(oo) = 1/2. Clearly, its response has also an exponential 
character. However, this happens to be the case of dipoles relaxing irrespec- 
tive of each other and of their environment. If the dipoles interact with their 
environment, and the interaction is complex (or random), their contribution 
in relaxation already will not result in any exponential delay. 

Assume that the interaction of dipoles with environment is taken into account 
with the aid of the temporal subordination. We will consider the evolution of 
the number of dipoles in the states 4- and 7. This are parent random processes 
in the sense of subordination. They may be subordinated by another random 
process with a probability density, say p(r, t). If n^{j) and n^(r), taking from 
Eq. (pQ) as probability laws of the parent processes, depend now on a local time 
r, then the resulting n^(t) and ni(i) after the subordination is determined by 
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the integral relation 



oo oo 



n t (t) = / n t (T)p(r,t)dr , n^t) = / n^r) p(r,t) dr 



o o 



If the directing process (a new time clock or stochastic time arrow [9]) is 
an inverse a-stable process, its probability density has the following Laplace 
image 



where Br denotes the Bromwich path. This probability density has a simple 
physical interpretation. It determines the probability to be at the internal time 
(or so-called operational time) r on the real (physical) time t. The function 
F a (z) can be expanded as a Taylor series. Besides, it has the Fox' H-function 
representation 



where T(x) is the ordinary gamma function. In the theory of anomalous dif- 
fusion the random process S(t) is applied for the subordination of Levy (or 
Gaussian) random processes [10]. The inverse a-stable process accounts for 
the amount of time that a walker does not participate in the motion process 
[TTj . If the walker participated all time in the motion process, the internal 
time and the physical (external) time would coincide. 

As was shown in [§11121113] . the stochastic time arrow can be applied to the 
general kinetic equation. Then the equation describing a two-state system 
takes the following form 



where D a is the a-order fractional derivative with respect to time. Here we 
use the Caputo derivative [H], namely 






< a < 1, 



(3) 




n — 1 < a < n, 
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where x^ n '(t) means the n- derivative of x{t). The relaxation function for the 
two-state system is written as 

0cc(*) = 1 " 2n ; (t) = 2n t (t) - 1 = E a (-2wt a ), 



where E a (z) = J2™=o z n /Y{l+na) is the one-parameter Mittag-LefHer function 
[T5] . It is important to notice that the relaxation function corresponds to the 
Cole-Cole (CC) law. With reference to the theory of subordination the CC law 
shows that the dipoles tend to equilibrium via motion alternating with stops 
so that the temporal intervals between them is random. The random values 
are governed by an inverse a-stable subordinator. 

The evolution of n^(t) and in Eq. (TSJ) can be connected with the Mittag- 
LefHer distribution. If Z n denotes the sum of n independent random values 
with the Mittag-LefHer distribution, then the Laplace transform of n- l / a Z n is 
(1 + s a /n)~ n , which tends to e~ sa as n tends to infinity. Following the argu- 
ments of Pillai [TB] , this supports an infinity divisibility of the Mittag-LefHer 
distribution. Due to the power asymptotic form (long tail) the distribution 
with parameter a is attracted to the stable distribution with exponent a, 
< a < 1. The property of the Mittag-LefHer distribution permits one to 
determine a stochastic process. The process (called Mittag-Leffler's) arises of 
subordinating a stable process by a directing (generalized) gamma process 
|. In this case the relaxation function has the form 

M^-t J^^ w^r^, (4) 



where a, b, thn axe constant. The one-side Fourier transformation of the re- 
laxation function gives the Havriliak-Negami (HN) law 

I \ 1 ~i"t ( ^ H n(*)\ 1 

XHNM = i e {-^i^) dt= (i + ^ry (5) 



This result also corresponds to the well-know HN empirical law. Thus, the HN 
relaxation can be explained from the subordination approach, if the hitting 
time process of dipole orientations transforms into the Mittag-LefHer process 
[T7] . For that the hitting time process has an appropriate distribution at- 
tracted to the stable distribution. The subordination of the latter results just 
in the Mittag-LefHer process. It is interesting to observe that the Levy pro- 
cess subordinated by another Levy one leads again to the Levy process, but 
with other index |18] . Unfortunately, the description from the Mittag-LefHer 
process gives nothing for the derivation of any macroscopic response equation 
like Eq.©. 
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3 Tempered a-stable Process and Its Inverse 

The relaxation model based on the inverse a-stable process starts with the 
consideration of a-stable processes having the infinite-moment difficult. To 
overcome it, one can develop an approach stated on tempered a-stable pro- 
cesses. The tempered a-stable process [191120] has the Laplace image of its 
distribution in the form 

f(u) = exp (S a - (u + 5) a ) . (6) 



When 5 equals to zero, the tempered a-stable process becomes simply a-stable. 

However, the distribution (J2]) describes only probabilistic properties in terms 
of internal time. For subordination we need the probability distribution of the 
inverse tempered a-stable process. If /(r, t) is the p.d.f. of internal time, then 
the p.d.f. of its inverse g(r,t) can be represented as 

9(r,t) = -— J f(t',r)dt'. 



Taking the Laplace transform of g{r,t) with respect to t, we get 

~9(r, u) = -~K*> r) = + e-^-^ . (7) 

UOT U 

When 5 -» 0, Eq. (JZJ) tends to 
~g(r,u)=u a ~ 1 e- Tua . 

This expression corresponds to the Laplace image of an inverse a-stable p.d.f. 
describing a directing process in the theory of Cole-Cole relaxation. After the 
inverse Laplace transform we have Eq.(Bl). 

In the Laplace space the function g(r,u) has a simple form. Because of gen- 
eral properties of the Laplace transform we can find g(r,t) explicitly, but its 
representation is enough complicated and expressed through a integral of the 
Wright functions [21] . For our calculation it will be sufficient to know only the 
function g(r,u). Therefore, we will not present any explicit form g(r,t) here. 
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4 Macroscopic Response Equation of Tempered Relaxation 

If in Eq.([T]) the value w will depend on time as aA a t a ~ 1 , we come to the 
description of the Kohlrausch- Williams- Watts (KWW) relaxation. Although 
such a equation does not contain any (for example, micro/mesoscopic and so 
on) details about relaxation processes, it is convenient for practical purpose 
because of its simplicity. When the relaxation follows the CC, CD (Cole- 
Davidson), HN laws, the equation ([T]) is not so simple as in the case of D 
and KWW relaxation. Recall that the CC relaxation and response functions 
can be expressed in terms of a solution of the fractional differential equation 
|18j . With macroscopic equations for the CD and HN responses the situation 
becomes more else complicated. Consider the derivation of the macroscopic 
response equation of tempered relaxation in more details. 

For a general type of a Markovian process the general kinetic equation is 



where W denotes the transition rate operator (see details, in |22j). This equa- 
tion defines the probability p for the system transition from one state into 
others. Next, we determine a new process governed by an inverse tempered 
a-stable process with the Laplace image (jTJ), namely 



dyit) 

dt 



Wp(t) 




oo 









The Laplace transform p a {s) is given by 



oo 








and leads to 




s 



(s + 5) a - 5' 



{[(s + 5r-6 a }p((s + 6r-6 a )-p(0)} 



s 
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From this it follows 

t 

p a (t)=p(0) + JdrM(t-r)W Pa (t), (10) 

Q 

where the kernel M(t) is written as 
M{t) = e~ 5t t a ~ l E a>a {5 a t a ) ) 

where E aj p(z) = J2™=q z n /T((3 + na) is the two-parameter Mittag-Lemer func- 
tion [15J . This equation covers a number of particular cases known earlier. For 
a = 1 we obtain Eq. (JSJ), and for 5 = it becomes fractional. If a system has 
discrete states, then the generating function is of the form 

oo 
k=0 



where ( takes values | £ |< 1 for a series to converge. With the help of the 
generating function, one can find the moments by taking the derivative with 
respect to ( and then setting £ = 1. The generating function of the process 
governed by the stochastic time clock is given by the relation 

oo 

Ga(C,*) = J g{r,t)G{C,r)dr. (11) 
o 

Thus, the generating function for a discrete Markovian process directed by a 
subordinated process can be obtained from the appropriate generating func- 
tion of the parent process by immediate integration. 

The relaxation in a two-state system under the inverse tempered a-stable 
subordinator gives 

t 

n t (t) = n t (0) + w J M(t - r) {n±(t) - n t (t)} dr, 
o 

t 

n l (t)=n i (0) + w J M(t-T){n^t)- H (t)}dr. (12) 
o 

For ( < 1 (or 8 — > 0) the kernel M(t) takes the power form t a /T(a) as a 
fractional kernel in the integral representation of Eq. ([3]). However, for t > 1 
(or a — > 1) M(t) becomes constant and, as a result, Eq. (1T2T) transforms into 
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Fig. 1. The relaxation of the part of dipoles in the state "f and the part n± of 
dipoles in the state \. (for w/S = 1/4, a = 0.75). 

the integral form of the ordinary equations (CQ). From the linearity of these 
equations it just follows 

nt(t) + n±(t) = nt(0) + 714.(0) , n t (t) - n ; (t) = 
t 

nt(0) - n^O) - 2wJ M(t - r) {n t (r) - n x (r)} dr. 


Consequently, we obtain 

t 

n t {t) = l-wj e~ 5r r^ 1 E a>a ([8 a - 2w}r a ) dr, 


t 

H (t)=w J e~ 5r r*" 1 E atCe {[6 a - 2w}r a ) dr. 


The equations have steady states (t — > 00) corresponding to equilibrium in 
this system. According to [23], we know 

00 

J e~ ax x a - 1 E a , a (±bx) dx = — — -, (Re(o) > |6| 1/q ), 
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then n^(oo) = n^(oo) = 1/2 for any 5 > and w > (see, for example, 
Fig. [1]). The transition rate w is again defined by microscopic properties of the 
system. Thus the relaxation function takes the form 

0tem P (t) = 1 - 2 H {t) = 1 - 2w J e- ST r^ 1 E a , a ([5 a - 2w}r a ) dr. 



The response function f te mp(t) = — d(j> temp (t) / dt is written as 
f temp (t) = 2we' St t a ~ 1 E a , a ([8 a - 2w]t a ) . 

Fig. [2] demonstrates how to change the response function with the increase of 
5. 

For the experimental study it is interesting to get the frequency- domain rep- 
resentation of the latter function. Its real part describes a dispersion of the 
relaxing medium, and its imaginary part is an absorption. The values explic- 
itly can be measured in experiments. The one-side Fourier transformation of 
the response function gives 

XtcmpM = Je-™ f- ^S. ) dt = -L r- , (13) 

A pV ; J \ dt J 1 - a a + (iw/wp + a) a K ' 

where co p = (2w) 1//a and < a = 5/{2w) l / a < oo are constant. The parameter 
Up is the characteristic frequency of the relaxing system. It is easy to notice 
that the expression (fT3"j) for a = 1 is reduced to the D law, for 5 = (or o = 0) 
it describes the CC relaxation, and for o = 1 it does the CD law (see Fig. [3]). 



5 Major Properties of Tempered Relaxation 

The frequency dependence of the dielectric susceptibility for orientational po- 
larization of dipoles has been the subject of experimental and theoretical stud- 
ies for many years, but still there is no any generally accepted theory capable 
of explaining the observed phenomena. Experimentally it is well known that 
the complex dielectric susceptibility x( w ) — x'i -*) — i~x"{u)) of most dipolar 
substances demonstrates a peak in the loss component x"(uo) at a characteris- 
tic frequency u p , and it is characterized by high- (cj ^> u p ) and low- frequency 
(to Up) dependencies. The tempered relaxation shows 
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A + Bcos(C) 
A 2 + 2AB cos(C) + B 2 ' 

sin(C) 
A 2 + 2AB cos(C) + 5 2 ' 

where A = 1 — cr a , 5 = (er 2 + cj 2 /uj 2 ) a ^ 2 and C = a arctan(u;/<5). For small 
a; and any positive a 7^ it is easy to see that lim^o Xtemp ( u ) ~ u an d 
lim^o Xte m p (w) ~ 1 whereas for large w we get lim^^ Xtemp M ~ and 
lim^oo x^mp M ~ Tllis implies that 

,. Xtemp M /0>7V\ , / 7T \ 

lim — — i-— - — = tan — = cot n— , 

where n — 1 — a, that is in agreement with the experimental results [T1I2] . 
However, for small u; we come to 

v Xtemp ( W ) 

lim — t— r — ; -; — - = OO . 

w - >0 Xtcmp(O) - Xtemp M 



Xtemp M = 
Xtemp M = 



This means that the energy lost per cycle does not have a constant relationship 
to the extra energy that can be stored by a static field. According to such an 
asymptotic behavior and in Fig. |3]it is seen that the tempered relaxation takes 
an intermediate place between the D, CC and CD types of relaxation. 



6 Conclusions 



In this paper we have represented our progress in the subordination analysis 
of relaxation phenomena in the complex systems. The approach permits ones 
to consider many relaxation laws on the unique theoretical base originating 
from the stochastic nature of relaxation. The general probabilistic formalism 
treats the relaxation of the complex systems regardless of the precise nature 
of local interactions. Following this approach, we have derived the empirical 
relaxation laws and their macroscopic equations, have characterized their pa- 
rameters, have connected the parameters with local random characteristics 
of the relaxation processes, have demonstrated how to make the transition 
from the microscopic random dynamics in the complex stochastic systems to 
the macroscopic deterministic description by integro-differential equations. It 
should be pointed out that the form of these equations is a direct sign of com- 
plexity evolution in such systems. Although we restricted only by the detailed 
analysis of two-state systems, due to Eq. flTU]) . this consideration can be devel- 
oped to the study of many-state systems (as an example, see the analysis of 
the three-state fractional system in [IB])). The tempered relaxation establishes 
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a connection between several types of relaxation (D, CC and CD). Its asymp- 
totic behavior earnestly shows that starting as a non-exponential relaxation, 
latter it tends to the D law. Moreover, the theory of subordination suggests 
a clear interpretation of the tempered relaxation. As applied to the dielectric 
relaxation, the interaction of dipoles with each other and their environment 
has a confined time of action on the relaxation process near a starting point 
of relaxation. Latter they behave independently just as this is the case of ex- 
ponential (D) relaxation. To put it in another way, the dipoles are strongly 
connected with each other in the initial stage of relaxation (and, therefore, 
their response function obeys a non-exponential decay), but the connection is 
not long-lived, and in an interval of time each dipole evaluates on its own. 
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